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D ' Abstract. Multiple initial state parton interactions in p(d)+A collisions are 

calculated in a Glauber-Eikonal formalism, which incorporates the competing pattern 



I of \ow-px suppression due geometrical shadowing, and a moderate-pT Cronin 

' enhancement of hadron spectra. Dynamical shadowing effects, which are not included 

in the computation, may be extracted by comparing experimental data to the baseline 
^ ' provided by the Glauber-Eikonal model. Data for tt^ production at midrapidity show 

absence of dynamical shadowing in the RHIC energy range, ^/s ~ 20—200 GeV. Recent 
preliminary data at forward rapidity are addressed, and their interpretation discussed. 
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. In proton (p), or deuteron (d), reactions involving heavy nuclei (y4~200) at y/s < 40 

^ ' AGeV, the moderate transverse momentum (pj'~2-6 GeV) spectra are enhanced rela- 
tive to linear extrapolation from p + p reactions. This "Cronin effect" P is generally 
attributed to multiple scattering of projectile partons propagating through the target 
^ ! nucleus [2.. In this talk, we discuss multiple parton scattering in the Glauber-Eikonal 
^ ' (GE) approach [Slllj, in which sequential multiple partonic collisions are computed in 

pQCD, and unitarity is naturally preserved. The low-p^ spectra in p + A collisions are 
suppressed by unitarity. At moderate pt, the accumulation of transverse momentum 
', leads to an enhancement of transverse spectra. At high pt the binary scaled p + p 
spectrum is recovered: no high-p^ shadowing is predicted in this approach. 

Hadron production in p-|-p collisions. The first step to understand p + A 
collisions is to understand p + p collisions. The pQCD formula for the single inclusive 
hadron transverse spectrum is: 
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dp^dy 



Here we considered only elastic parton-parton subprocesses, which contribute to more 
than 98% of the cross section at midrapidity. In Eq. ((T]), 

{xfi/p)y.,PT = -^^ il^.. j dy2Xifi/p{xi,Ql)^ {s,i, u) X2fj/p'{x2,Ql) J^^^y, (2) 
=^^TT6-J dy2^is,t,u)x2fM^2,Q,) (3) 
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are interpreted, respectively, as the average flux of incoming partons of flavour i from 
the hadron p, and the cross section for the parton-hadron scattering. The rapidities of 
the i and j partons in the final state are labelled by Ui and y2- Infrared regularization is 
performed by adding a small mass to the gluon propagator and defining the Mandelstam 
variables i(u) = — mf.(l + 6"^^'^^^), with = \/pt + Po- -^^^ more details, see Ref. 
Finally, inclusive hadron production is computed as a convolution of Eq. (|T)) with a 
fragmentation function Di^h{z,Ql^). 

In Eqs. (0)-©, we have two free parameters, po and the K-factor K, and a some- 
what arbitrary choice of the factorization and fragmentation scales, Qp = Qh = mx/2. 
After making this choice, we fit po and K to hadron production data in pp collisions at 
the energy and rapidtiy of interest. Equation ((T)) is very satisfactory for Qt ^ 5 GeV, but 
overpredicts the curvature of the hadron spectrum in the gT=l-5 GeV range. This can 
be corrected for by considering an intrinsic transverse momentum, kx, for the colliding 
partons [Sj. We found that a fixed (k^) = 0.52 GeV leads to a dramatic improvement 
in the computation, which now agrees with data at the ±40% level [Sj. Finally, we 
obtain, at midrapidity 1] = 0, po = 0.7 ± 0.1 GeV and K = 1.07 ± 0.02 at Fermilab, and 
Po = 1.0 ± 0.1 GeV and K = 0.99 ± 0.03 at RHIC. 

Prom p+p to p+A collisions. Having fixed all parameters in p+p collisions, 
and defined the parton-nucleon cross section (jH)), the GE expression for a parton-nucleus 
scattering is 0: 
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x:|^r,(i,)x...x^r.(6)e--"<-)--") , (4) 

where 7^(6) is the target nucleus thickness function at impact parameter b. The 
exponential factor in Eq. ^ represents the probability that the parton suffered no 
semihard scatterings after the n-th one, and explicitly implements unitarity at the 
nuclear level. Assuming that the partons from A suffer only one scattering on p or 
d, we may generalize Eq. as follows, without introducing further parameters: 

d^PTdyd^h 1 "^^^^'p'y^'P^ d^pTdyd^h ytv ^ V ^"^^^Z^^^^'^" V% 



Unitarity introduces a suppression of parton yields compared to the binary scaled 
p+p case. This is best seen in integrated parton yields: da"^^ / dyd% ~ 1 — e "'^'^'^^"^"^-^'^^^ 
At low opacity x = (^^^ {Po)TA{b) <C 1, i.e., when the number of scatterings per parton 
is small, the binary scaling is recovered. However, at large opacity, x ^ 1; the parton 
yield is suppressed: da^^ / dyd% ^ 1 < cr*^(po)^A(&)- This suppression is what we 
call "geometrical shadowing" , since it is driven purely by the geometry of the collision 
through the thickness function Ta- As the integrated yield is dominated by small 
momentum partons, geometrical shadowing is dominant at low px- 
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Figure 1. Cronin effect on pion production at Fermilab and RHIC ^Hl at 77 = 0. 
The solid curve is the GE computation. Theoretical errors due to the fit of Pq are 
shown as a shaded band around the solid curve. The rightmost panel shows the 0- 
20%/60-88% centrality classes ratio. 

Beside the geometrical quark and gluon shadowing, which is automatically included 
in GE models, at low enough x one expects genuine dynamical shadowing due to non- 
linear gluon interactions as described in Colour Glass Condensate (CGC) models [Oj. 
However, it is difficult to disentangle these two sources of suppression of spectra, and 
to understand where dynamical effects begin to play a role beside the ubiquitous geo- 
metrical effects. The GE model computation outlined above can be used as a baseline 
to extract the magnitude of dynamical effects by comparison with experimental data. 



Cronin effect at Fermilab and RHIC. The Cronin effect may be quantified 
by taking the ratio of hadron px spectrum in p(d)+A collision, and dividing it by the 
binary scaled p+p spectrum: 

where b is the average impact parameter in the centrality bin in which experimental 
data are collected (see 0). The GE model reproduces quite well both Fermilab and 
PHENIX data at = (Fig.l, left and center). It also describe the increase of the Cronin 
effect with increasing centrality (Fig.l, right). If dynamical shadowing as predicted in 
CGC models was operating in this rapidity and energy range, the central/peripheral 
ratio should be smaller than the GE result: the more central the collision, the higher 
the parton density in the nucleus, the larger the non-linear effects. Therefore, we 
conclude that there is no dynamical shadowing nor Colour Glass Condensate at RHIC 
midrapidity. 

To address the preliminary BRAHMS data at forward rapidity 77 ~ 3.2 [11], 
we would ffist need to fit po and K in p+p collisions at the same pseudo-rapidity. 
Unfortunately the available data px-Tange 0.5 ^ Pt ^ 3.5 Ge is not large enough px 
for the fit to be done. Therefore, we use the parameters extracted at rj = 0. The 
resulting Cronin ratio, shown by the solid line in Figure 2, overestimates the data at 
such low-px- This may be corrected in part by considering elastic energy loss [T^ . 
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Figure 2. Cronin ratio at rj — 3.2. 



Furthermore, the opacity Xo = 0.95 might be 
underestimated, due to the use of the mid-rapidity 
parameters: to check this we tripled the opacity, 
but the resuhing dashed hne still overestimates the 
data. 

Does the discrepancy between our calculations 
and the BRAHMS data prove that a CGC has been 
observed? It is too early to tell. An important 
observation ^Ij is that the yield of positive 
charged hadrons exceeds the yield of negative 
charged hadrons at moderate pt by just the factor 
expected in HIJING due to dominance of valence 
quark fragmentation effects: gluon fragmentation 
is subdominant mechanism in this kinematic range. In addition, current simplified CGC 
scenrarios predict asymptotic l/p^ absolute behavior which overestimates by an order 
of magnitude the observed absolute cross sections. In our pQCD approach the absolute 
cross sections in both p+p and d+Au collisions are much closer to the data, though the 
modest discrepancies shown in Figure 2 remain. Gluon shadowing appears to be needed 
in this X ~ 10~^ — 10^^ regime, though more data will be required to quantify the effect. 

The optimal region to study the onset of dynamical shadowing is 3 < ^ 6 GeV, 
where the GE model expects the peak of the Cronin ratio to be. If in this region the final 
data at forward rapidity < 77 < 4 explored by the four RHIC collaboration will reveal 
a consistent pattern of suppression compared to the GE computation (as preliminary 
data seem to suggest), and nonperturbative fragmentation effects will not be able to 
explain it, then the case for the CGC will be made more solid. The case would be even 
stronger if at the same time a progressive disappearance of back-to-back jets in favour 
of an increased dominance of monojet production was observed. 
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